Creation of discrete solitons and observation of the Peierls-Nabarro barrier in 

Bose-Einstein Condensates. 
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We analyze the generation and mobility of discrete solitons in Bose-Einstein condensates confined 
in an optical lattice under realistic experimental conditions. We discuss first the creation of fD 
discrete solitons, for both attractive and repulsive interatomic interactions. We then address the 
issue of their mobility, focusing our attention on the conditions for the experimental observability 
of the Peierls-Nabarro barrier. Finally we report on the generation of self-trapped structures in 
two and three dimensions. Discrete solitons may open alternative routes for the manipulation and 
transport of Bose-Einstein condensates. 



I. INTRODUCTION 

The experimental achievement of the Bose-Einstein 
condensate (BEC) Q has outbursted an extraordinary 
interest within the last years in the physics of ultracold 
atomic gases. This interest can be partially explained by 
the inherently nonlinear character of the BEC physics in- 
duced by the interatomic interactions. At sufficiently low 
temperatures, the physics of the condensates is governed 
by a nonlinear Schodinger equation with cubic nonlinear- 
ity, also called Gross-Pitaevskii equation (GPE), similar 
as that encountered in other physical systems, as e.g. 
nonlinear optics (NLO) in Kerr media. The analysis of 
the resemblances between BEC physics and NLO has lead 
to the rapidly-developing field of Nonlinear Atom Optics 
(NLAO) 0- Recently several experiments have high- 
lighted various NLAO phenomena, as dark solitons in 
BEC with repulsive interatomic interactions bright 
solitons in ID BECs with attractive interactions H,l1> 
and condensate collapse 0. 

During the last few years, the possibility of loading 
a BEC in an optical lattice formed by a laser standing 
wave has attracted considerable attention, mostly mo- 
tivated by the close resemblance between these systems 
and solid-state devices. In this sense, several remarkable 
experiments have been recently reported, as the observa- 
tion of Bloch oscillations of BECs |8|, |9| , the realization 
of Josephson-junction arrays of BECs placed in different 
lattice sites or even the achievement of the super- 
fluid to Mott-Insulator transition Recently several 
nonlinear BEC phenomena have been analyzed in the 
presence of optical lattices, as the dynamical superfluid 
to insulator transition [12], the BEC transport in the 
presence of dispersion managing [p3 | , and the generation 
of gap solitons, i.e. bright solitons with condensates with 
repulsive interactions 0] . 

Particular interest has been recently devoted to those 
phenomena occurring when the condensate dimensions 
become comparable to the lattice wavelength. In that 
case, the discrete structure introduced by the lattice 
potential may lead to similar phenomena as those ob- 
served in NLO in periodic structures. In particular the 
analysis of discrete solitons (DS) 0, UtI Hsj in the 



BEC context has recently attracted a growing attention 
0> HE HH HE HH • Specially interesting phenomena ex- 
clusively induced by the discreteness of the system have 
been analyzed in the NLO context, such as the restriction 
of the mobility of the DSs due to the so-called Peierls- 
Nabarro (PN) barrier 1241 . or the possibility to generate 
two-dimensional DSs 1251. 



Although several properties related with DSs in BEC 
have been already reported 0, HE , the realization 
of DSs under realistic experimental conditions has so far 
not been analyzed in detail. Therefore, one of the aims of 
the present paper is to discuss the creation of these struc- 
tures in the frame of the recent experiments on BEC in 
optical lattices. In particular, we shall discuss the gen- 
eration of ID DS for both attractive and repulsive in- 
teracting condensates. Once created, the effects of the 
discrete nature of the DSs in BEC should be analyzed 
by means of the observation of the PN barrier for its 
mobility. A second aim of the present paper is to dis- 
cuss the conditions for the experimental observability of 
this barrier. Interestingly, the PN barrier is largely over- 
estimated within the usual tight-binding approximation 
even under conditions for which this approximation is 
typically assumed. Finally, in the last part of our pa- 
per, we discuss the possibility of achieving 2D and even 
3D self-trapped structures, which could offer alternative 
routes for the controllable manipulation of BECs. 



The scheme of the paper is as follows. In Sec.[n]we dis- 
cuss the physical system under consideration as well as 
the basic equations to describe it. Sec. lIIII presents a vari- 
ational approach which allows for an analysis of discrete 
structures in arbitrary dimensions. Sec. IIVI is devoted to 
the analysis of the generation of ID DSs, for both attrac- 
tive and repulsive interatomic interactions. We address 
also in Sec. IIVI the issue of their mobility and provide 
the conditions for the observability of the Peierls-Nabarro 
barrier. In Sec. we analyze the creation of 2D and 3D 
self-trapped structures. We finalize in Sec. IVII with our 
conclusions. 
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II. PHYSICAL SYSTEM 

In the following we consider a trapped BEC in the 
presence of an optical lattice. Th e p eriodic structure 
leads to an energy band structure [2(J |2?], IH |2i| ^ind 
strongly modifies the dynamics of the condensate H H 
El El Hi III El HI In the mean field approximation, 
the full BEC dynamics (at a temperature much smaller 
than the critical one for the condensation) is governed by 
the time-dependent GPE: 

ih <m^A = (-^—A + V (r) + «# (r, t) \ 2 ) iP (r, t) , 
at \ 2m J 

(1) 

where g = Airh 2 a/m, with a the s-wave scattering length 
and m the atomic mass. The condensate wave function 
is normalized to the total number of particles N . The 
external potential is given by: 

7(f) = |(o; x V+a; J/ V+a; z V) + 
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(2) 



which describes both the magnetic trap potential and the 
optical lattice (created by two counter propagating laser 
beams of wavelength A along each axis). The angular 
frequencies of the magnetic trap in each direction are 
denoted by u)j . The optical lattice is characterized by its 
depth Vb and by its lattice period d — A/2, which defines 
the so-called recoil momentum k r = ^. In the following, 
and following the standard notation, we refer the depth 
of the optical potential in units of the so-called recoil 
energy E r = h 2 k 2 /2m. 

When the optical depth of the lattice is much larger 
than the chemical potential (Vb >> fi), and the system 
can be considered as confined within the lowest energy 
band, one can employ the tight-binding approximation 
and rewrite the condensate order parameter as a sum 
of wavefunctions localized in each well of the periodic 
potential: 



4,(r,t) = VNJ2^ (*) <Pn (f) 



(3) 



where tp n (r) = ip (r — r n ) denotes the on-site wavefunc- 
tion. By inserting the Ansatz J3J into Eq. (JJJ , one obtains 
that the GPE reduces to a discrete nonlinear Schrodinger 
equation (DNLSE): 

The dynamics of the system depends mostly on the in- 
terplay between the tunneling rate (J) and the two body 
interactions (non linear energy, U). The tunneling rate 
can be expressed as: 



h 2 - _ 

J = - j dr[—S7p n Vip n+ i + Lp n V(r)Lp n+1 ] 



the nonlinear term acquires the form: 
U = gN [ dr^ n 



(6) 



and the on-site energies are given by: 

en = J dr[^Vn) 2 + V{r) V l\. (7) 

In order to calculate the value of the nonlinear energy 
Eq. (BJ, we use a Gaussian Ansatz for the wavefunction 
on-site, where the width is obtained by minimization of 
the energy |31j . In the pure one-dimensional case, we use 
the interaction constant obtained by averaging the 3D 
coupling constant over the radial density profile . To 
calculate the tunneling rate, the same Gaussian Ansatz 
can be employed in Eq. (JSJ or, in the limit Vb >> E r , 
it can be obtained from the exact result for the width of 
the lowest band in the ID Mathieu-equation |33| . 



III. VARIATIONAL CALCULATION 

Discrete solitons are characterized by being stable so- 
lutions of the Hamiltonian which propagate without dis- 
tortion. They correspond to minima of the energy of 
the system which in the tight-binding approximation ac- 
quires the form: 

oo 

E = E {-JK^n-i+K + i) + e n \cj )n \ 2 + -\^ n \ 4 }. 
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(8) 

Let us consider first the ID case. To find the minima of 
the energy, an exponential Ansatz for the soliton envelope 
given by <j> n = C exp (— j3\n\) can be employed, where C 
is a normalization constant and (3 is a variational param- 
eter which accounts for the inverse width of the soliton. 
Introducing this Ansatz in Eq. (JSJl the expression for the 
energy, in the negative scattering length case, becomes: 
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(5) 



By minimizing E with respect to the inverse width j3 one 
obtains the energy of the discrete structure centered in 
one minimum of the optical lattice (Fig^a)). 

Notice that in order to ensure that the discrete struc- 
ture corresponding to this minimum is indeed a soliton, 
one has to address the issue of its mobility. In contrast 
with the case of continuous solitons where applying an 
external momentum results in a linear response of the 
soliton, in discrete systems a similar effect occurs only 
for broad soliton distributions (i.e., those that occupy 
many sites). Conversely if the dimensions of the soliton 
are of the order of the lattice wavelength, A, the discrete- 
ness of the system begins to play a fundamental role. In 
particular, the discreteness generates an effective peri- 
odic potential energy, whose amplitude is the minimum 
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FIG. 1: Density profile of a discrete soliton (solid line) with 
(a) the center of mass centered in one minimum of the optical 
lattice and (b) the center of mass displaced by half lattice 
period. 



barrier which must be overcome to translate the center 
of mass of the system half a lattice period, i.e. from 
one minimum of the lattice (Fig^a)) to a neighboring 
lattice maximum (Fig^b)). This is the previously men- 
tioned PN barrier [24J . The energy of the state depicted 
in Fig^b) can be also calculated within a variational 
approach using <p n = C exp (— 0\n — 5I) as an Ansatz: 



E 

W\ 



2J(l + sinh(^)) l^-l) 



\U\ 



4 (e 2 / 3 + 1) ' 



(10) 



Again by minimization one obtains the energy of the dis- 
placed discrete structure. The difference between both 
energies (minima of Eq.Q and Eci. (|10|l ) corresponds to 
the PN barrier. 

Notice that the barrier becomes relevant when the soli- 
ton structure occupies few sites of the lattice. By increas- 
ing the number of occupied states the above two modes 
approach in energy and the barrier decreases approach- 
ing zero as the number of sites grows. In this sense, the 
PN barrier is a distinctive discrete phenomena. 

The above analysis can be straightforwardly general- 
ized to higher dimensions. Assuming the most general 
case, in which the width of the exponential Ansatz along 
the direction of the movement is different from the other 
directions, it can be shown that localized structures are 
also possible in two and three dimensions under an ap- 
propriate ratio between tunneling and nonlinear energy 
(see Sec. EJ- A variational calculation using a Gaussian 
Ansatz instead of the exponential one has also been per- 
formed obtaining equivalent results. 



IV. GENERATION AND MOVEMENT OF 
ONE-DIMENSIONAL DISCRETE SOLITONS 

In this section we analyze the issue of the generation 
and mobility of DSs in ID BEC under realistic exper- 
imental conditions. Let us first discuss the generation 
of DSs in condensates with positive scattering length. 
In particular, we consider a 87 Rb condensate hold in a 
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FIG. 2: (a)Energy of the first band in units of the recoil energy 
and (b)velocity profile in units of function of 

the quasimomentum in the first Brillouin zone for an optical 
potential depth of 8E r 



magnetic trap and in the presence of an optical lattice. 
As discussed in 0, |2j| it is possible to generate DSs in 
pure ID repulsive condensates when the tunneling rate 
balances the nonlinear energy of the system. For pos- 
itive scattering lengths, the compensation of these two 
effects is not possible unless the system has a negative 
effective mass m* (^ = -p-^r where E is the energy 
of the first band and k its quasimomentum (see Fig.EJ). 
At least two possible mechanisms can easily place the 
system in such a region, inverting thus the sign of the 
tunneling. One possibility jl-'"'J| consists on providing the 
condensate with an external momentum to place it at 
the edge of the first Brillouin zone where the effective 
mass is negative (see Fig. |SJa)). This can be done for 
instance, by introducing -in the absence of the magnetic 
trap - a tilt in the optical lattice. A second mechanism 
to reach the negative effective mass region relies on the 
variation of the relative phase of the condensate in the 
lattice. Concretely, a repulsive condensate in a periodic 
potential with a phase difference of n between consecutive 
wells gives rise to the so-called staggered discrete soliton 
type [2J|. Such a phase structure can be achieved us- 
ing the well established method of phase imprinting [3{| 
which allows to modify the phase of a condensate with- 
out modifying its density profile. To this aim we propose 
to use a second optical lattice with double spatial period 
than the first one acting for a time much shorter than 
the characteristic times of the system, i.e. the correla- 
tion (t c = ~) and the tunneling (r t = j) time. The 
phase imprinted in this way depends solely on the am- 
plitude of this second standing wave and on the time in 
which it acts. In the following we describe in detail how 
the phase imprinting method is implemented. Firstly 
we calculate the ground state of the system in the pres- 
ence of an optical lattice. This can be done either by 
calculating directly (using GPE in imaginary time) the 
ground state of the system in the presence of both the 
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FIG. 3: (a) Density profile of the ground state of a 87 Rb 
condensate in the combined (magnetic and optical) trap for 



N = 2000, w x = 50x2?r Hz, u y 



92x2tt Hz, a = 5.8nm, 



Vo = lE r and d — 0.8/im.(b) Density profile of the staggered 
soliton 100ms after the imprinting of a phase difference of 
7r between consecutive wells and after the magnetic trap is 
switched off.(c) Phase profile of the staggered soliton shown 
in (b) 



magnetic trap and the optical lattice or by calculating 
the ground state of the system in the presence of the 
magnetic trap only and afterwards growing adiabatically 
the optical lattice and letting the system to evolve to the 
new ground state. As expected both methods yield the 
same ground state. Once the ground state is found, a 
second optical lattice with an amplitude of 72.5E r act- 
ing for t = 0.4ms performs the phase imprinting while the 
magnetic trap is suddenly switched off. Our results are 
summarized in Fig. EI Fig-BI a ) shows the density profile 
of the ground state of the combined trap (magnetic and 
lattice) while Fig.|3Jb) displays the density profile of the 
localized structure 100ms after turning off the magnetic 
trap. The system evolves from the ground state shown 
in Fig. E^a) to the localized structure shown in Fig. |3£b) 
which remains unaltered for times much longer than the 
tunneling time of the system. This structure contains 
40% of the initial number of atoms. In Fig E^c), we 
display the phase profile corresponding to Fig. \3j[b) in 
where clear phase jumps of 7r between consecutive sites 
of the optical lattice in the spatial region occupied by 
the localized structure are present. Finally by applying 
an external momentum to the localized structure we ob- 
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FIG. 4: Momentum associated with the PN barrier calculated 
with the variational method (open dots), with the imaginary 
time evolution of the GPE (squares) and with the DNLSE 
(stars). The necessary momentum to move the soliton (cal- 
culated with the real time evolution of the GPE) is displayed 
by the dotted line. 



serve that it moves without distortion, evidencing thus 
that such structures correspond indeed to discrete bright 
solitons. Notice that the applied momentum must keep 
the structure in the negative effective mass region of the 
band (Fig.^a)) in order not to destroy the soliton. 

We turn now to the case of negative scattering length. 
In continuous systems, the attractive nature of the inter- 
actions compensates the effect of the kinetic energy and, 
therefore, the ground state of a one dimensional homo- 
geneous condensate with negative scattering length is al- 
ready a bright soliton. The presence of an optical lattice 
in such systems permits the creation of a discrete bright 
soliton. Thus, for a fixed number of atoms, N, one can 
straightforwardly vary the number of occupied sites by 
varying only the depth of the optical potential, Vo- Since 
the interactions are now attractive the number of sites oc- 
cupied by the discrete solitons is typically much smaller 
than in the case of positive scattering length. Bearing all 
this features in mind, it becomes clear that such systems 
are unique to study how the movement of a discrete soli- 
ton depends on the degree of discreteness of the structure, 
i.e., to observe the PN barrier |2^|. We have analyzed the 
case of a condensate of 7 Li confined in a magnetic trap 
with uj y — uj z = 20 x 2tt kHz, and u x — 75 x 27r Hz. We 
fix aN = — 14nm |3(j. Again, to analyze the generation 
and the propagation of discrete structures we start by 
calculating the ground state of the system in the com- 
bined (lattice + magnetic trap) potential. We proceed 
as in the case of Rubidium, either growing adiabatically 
the optical lattice after the ground state of the conden- 
sate in the trap has been found or by finding directly the 
ground state in the presence of both potentials. After 
the magnetic trap is switch off we observe that the local- 
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ized structure remains without dispersion for times much 
larger than the tunneling time of the system. To study 
the effect of the discreteness we vary the values of the 
depth of the optical potential between 6 < Vo/E r < 10 
for a fixed lattice spacing of d = 0.8/im. The number of 
sites significantly occupied by the discrete structure vary 
then between 11 and 3. 

We calculate for these sets of parameters the PN bar- 
rier, i.e., the difference in energy between the localized 
structure centered in a minimum of the periodic poten- 
tial (Fig.QIa)) and the structure corresponding to a lat- 
tice displaced by half lattice period (Fig. ^b)). To test 
the validity of the variational method and the validity 
of the tight binding approximation, we calculate the PN 
barrier using the expressions developed in Sec. IIIII and 
solving numerically in imaginary time the DNLSE (us- 
ing a Runge-Kutta method). Then, we compare these 
results with the barrier obtained by solving directly the 
full GPE in imaginary time. The results of those calcu- 
lations are summarized in Fig. where we display the 
momentum associated with the PN barrier (in units of 
k r ) as a function of the depth Vq of the optical lattice. 

We observe that for large optical potential depths 
(Vq > 8E r ) the results obtained with the DNLSE and the 
variational method using the exponential Ansatz match 
each other perfectly and agree well with those obtained 
by solving directly the GPE. For low optical potential 
depths (Vo < 8E r ) the DNLSE and the variational start 
to disagree because the exponential Ansatz becomes less 
appropriate as the number of sites increases. Moreover, 
in this region a clear disagreement between the GPE and 
the methods that assume tight binding appears, evidenc- 
ing thus the inapplicability of the tight binding approxi- 
mation. Note that, surprisingly, the tight-binding results 
for the PN barrier fail to match with the GPE solution 
even for values of the lattice potential that are typically 
assumed as being well into the tight-binding regime. In 
this sense our calculations reveal that the PN barrier 
is certainly very sensitive to slight deviations from the 
tight-binding conditions. The latter can be explained 
taking into account that if the on-site dynamics inside 
each potential well is not completely frozen, as assumed 
in Eq. ((3J), the energy associated to the on-site movement 
can smear out the PN barrier. 

Let us recall that physically the presence of the PN 
potential implies that the soliton will only move if a mo- 
mentum above the critical one, determined by the PN 
barrier, is provided. Therefore, we perform a dynamical 
study by solving the GPE in real time after providing 
an instantaneous transfer of momentum to the structure. 
The latter can be achieved either via phase imprinting or 
by applying a linear potential during a time shorter than 
any other time scale. We monitor for every value of Vq 
the minimum applied momentum to move the localized 
structure. The results are displayed in Fig. ||as a dot- 
ted line. We observe that the barrier calculated with the 
static GPE and the dynamical results only agree for low 
optical potential depths (Vq < 7E r ). For such cases the 



soliton is spread considerably in real space being, there- 
fore, well localized in momentum space. This is the nec- 
essary condition to move the structure as a whole in the 
linear region of the velocity profile (Fig. |2Ia)) where the 
velocity is proportional to the given momentum. The 
movement for the case of Vo = 6E r after a transfer of 
momentum of 0.1k r is depicted in Fig. EJa). Conversely, 
a clear deviation from all the previous calculations of 
the PN barrier appears as the lattice potential depth 
increases. Our results show that there exist regions in 
which it is not possible to reach dynamically the configu- 
ration corresponding to Fig. ^b) . This is what happens 
for large optical potentials (Vq > 8E r ) where the soliton 
is highly localized in space and, therefore, spread in mo- 
mentum. We observe that for (Vo > 10E r ), no motion 
is found no matter how large is the initial momentum 
given to the soliton. In this case, the width in momen- 
tum space of the localized structure is of the order of the 
momentum of the lattice, i.e., the first Brillouin zone is 
saturated and the movement is prevented. For interme- 
diate optical potentials (8 — 9E r ), we observe that move- 
ment of these structures is not possible in the linear part 
of the velocity profile (Fig. |2J) even if one overcomes the 
PN barrier. Nevertheless, if the given momentum is large 
enough, we observe in our simulations a reorganization 
of the structure which allows it to move by losing first a 
fraction of the atoms and spreading afterwards in space, 
being later on enough localized in momentum space to 
move as a whole (Fig. QJb)). The re-arranged structure 
exhibits a complex dynamics as a result of the interplay 
between nonlinearity and band structure. 



V. CREATION OF TWO- AND 
THREE-DIMENSIONAL SELF-TRAPPED 
STRUCTURES 

In this section, we briefly address the experimental pos- 
sibilities for the generation of discrete solitons in higher 
spatial dimensions in attractive condensates [22t l23| . 
From the variational Ansatz discussed in Sec. IIIII it is 
possible to obtain that two, and three-dimensional self- 
trapped structures may be stable provided that the ratio 
\J/U\ is low enough. The threshold value for the 2D 
case is 0.175(0.2) and for the 3D case 0.13(0.15) using an 
exponential (Gaussian) Ansatz. The limitation to such 
low tunneling rates imposes a restriction in the maximal 
number of sites occupied by the localized structure. 

We consider the full 3D case with an spherical mag- 
netic trap and we address the creation of 3D localized 
structures in the presence of a 2D or 3D optical lattice. 
Notice that due to the attractive interactions, the possi- 
bility of collapse imposes some restrictions in the number 
of atoms and in the features of the optical lattice. In some 
situations, even if the collapse in the magnetic trap is pre- 
vented, the adiabatic growing of the lattice leads to the 
on-site collapse. This effect must be taken into account if 
the 3D localized structures are to be created. A conden- 
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FIG. 5: (a)Density profile of a discrete soliton created in the 
presence of an optical potential of depth 6E r before and 21ms 
after a transfer of momentum of 0.1 fc r is applied to the struc- 
ture, (b) Density profile of the localized structure generated 
with a lattice of amplitude 8E r (dotted line) and the corre- 
sponding density profile 10ms after a kick of 0.3fc r is given 
(solid line). 77% of the atoms remain in the reconfigurated 
structure. 



sate of 7 Li with aN — — 70nm and initially in a spherical 
magnetic trap of frequencies (jJ x — {jjy — (jj z = 375 x 2vrHz 
has been considered. This situation corresponds for in- 
stance to N = 500 7 Li atoms and a modified s-wave 
scattering length a = — 0.14nm available with Feshbach 
resonances [gj^. Fig. [^displays the spatial structure of 
a matter discrete localized structure created in the pres- 
ence of a 3D optical lattice with period d = 1.6/xm, and 
V = 6E r . 

Therefore, the discreteness of the lattice potential al- 
lows for the interesting possibility of a controllable gen- 
eration of a self-trapped regular 3D BEC structure. Un- 
fortunately, the restriction in the number of sites occu- 
pied in the localized structures created in the two and 
three dimensional cases, imposes a severe limitation for 
their mobility. As already discussed for the ID case, the 
spread in momentum space of these structure prevents 
their motion, at least in the linear region of the band 
structure. 



FIG. 6: Spatial structure of a matter discrete localized struc- 
ture created in the presence of a 3D optical lattice with a 
depth Vb = 6E r , and a period d = 1.6/xm, as a function of 
two spatial coordenates, (x,y), for a fixed value of the third 
coordinate, z=0. 



VI. CONCLUSIONS 

We have analyzed the conditions to generate discrete 
solitons in ID Bose-Einstein condensates, either with 
positive or negative scattering length. In particular, in 
repulsive interacting condensates, the phase imprinting 
method has been proposed to controllably create bright 
staggered type solitons. Once generated, we have ad- 
dressed the mobility of these structures. This mobility 
is characterized by two different effects: (i) the presence 
of the PN barrier, a purely discrete effect which sets a 
minimal kinetic energy to move half a lattice period, and 
(ii) the spreading of the atomic wavefunction in momen- 
tum space due to the spatial localization in few lattice 
sites. The mobility is only possible if the PN barrier is 
overcome, but even if this is the case a clear soliton mo- 
tion is only possible if the system is placed in a region 
of linear dispersion. Our analysis shows that the esti- 
mation of the PN barrier is crucially sensitive to slight 
deviations from the tight-binding conditions. In partic- 
ular, the tight- binding approximation has been shown 
to fail significantly for large lattice potentials which are 
typically assumed to guarantee the validity of such an 
approximation. 

The mobility of discrete solitons generated in this way 
offer interesting possibilities in the context of BEC guid- 
ing. In this sense similar ideas as those analyzed in the 
case of optical DS [37j could be employed to generate DS 
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networks, which could consitute a novel (non dispersive) 
approach to the issue of integrated atom optics. 

Finally, we have shown that the discreteness of the lat- 
tice also allows for self-trapped structures in 2D and 3D. 
Due to the strong localization required, their mobility 
is prevented to a large extent. These higher-dimensional 
structures could open alternative routes to optical tweez- 
ers [l^] and magnetic conveyor belts [3jJ , for the storing, 
manipulation and transport of atomic condensates. 
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